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ABSTRACT 

We investigate the evolution of small perturbations around charged black 
strings and branes which are solutions of low energy string theory. We give 
the details of the analysis for the uncharged case which was summarized in a 
previous paper. We extend the analysis to the small charge case and give also 
an analysis for the generic case, following the behavior of unstable modes as the 
charge is modified. We study specifically a magnetically charged black 6-brane, 
but show how the instability is generic, and that charge does not in general 
stabilise black strings and p-branes. 
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1) Introduction 

Black holes have always been a source of fascination and debate, from the nature of 
the event horizon to the nature of the singularity and whether a full theory of quantum 
gravity can avert the apparent loss of unitarity in black hole evaporation. Classically, four 
dimensional black holes are stable, once formed they settle down to a state described solely 
by their mass, charge and angular momentum; therefore the singularities remain hidden 
from distant observers. Quantum mechanically, black holes are quite different objects, 
analogous to a thermal system. They have an entropy (proportional to surface area) 
and a temperature, since like black bodies they radiate thermally 1 . However, Hawking 
conjectured 2 that a black hole formed from a pure quantum state would radiate away 
leaving a mixed state of radiation; this would violate quantum mechanical unitarity. It 
is certainly true that Hawking's semi-classical description will break down near the final 
stage of black hole evaporation as planckian curvatures are approached, but if quantum 
gravity preserves unitarity and information is to be returned, it must do so well before this 
stage otherwise there is simply not enough energy left in a planck mass black hole to emit 
all the information stored in a macroscopic black hole. Therefore, in order to resolve this 
tension between quantum mechanics and thermodynamics, it seems that a semi-classical 
analysis of sub-Planck size black holes is needed. 

Most recently, attention has focussed on low energy string gravity and its implications 
for black holes. Some of these developments have been quite interesting. In Einstein grav- 
ity, charged black holes (the Reissner- Nordstrom solutions), in addition to having an outer 
event horizon, have an inner Cauchy horizon which is unstable to matter perturbations in 
the exterior spacetime 3 . However, there is no static charged black hole solution in Einstein 
gravity with only one horizon and a spacelike singularity. On the other hand, a key feature 
of low energy string gravity is the presence of a dilaton, which greatly changes the causal 
structure of charged black holes making them like Schwarzschild with one event horizon 
and a spacelike singularity 4 . This structure is generic, even if the dilaton has a mass 5 , as 
it must to keep in line with the principle of equivalence. 

Of course, most of the analysis of stringy black holes has been performed in low dimen- 
sions, namely two or four, whereas string theory tells us there should be ten dimensions, 
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ideally therefore, one should be examining black holes in ten dimensions. There has been 
work on black holes in higher dimensions 6 for a range of horizon topologies 7 . In four di- 
mensions, an event horizon must be topologically spherical 8 , but in higher dimensions this 
is not necessarily the case, we could have S' 2 xIR 6 , or 5' 3 xIR 5 topologies for the horizon. 
In a previous letter we pointed out that a large class of these black holes are unstable, 
namely the uncharged ones. In this paper, we present the details of this original argument, 
as well as providing an extension to cover charged black p-branes. 

Why should black holes be stable, yet black strings, say, unstable? Before answering 
this question in detail, it is worth examining a couple of naive arguments. Without loss 
of generality, consider a five dimensional black string in Einstein gravity, SCI14XIR. Then 
the equation governing the metric perturbation 5g ab = h ab , the Lichnerowicz equation, is 
essentially a wave equation 

A L h ab = (S c jtn + 2R a \ d )h cd = 0. (1.1) 

Because of the symmetries of the background Sch 4 xlR metric, this reduces to a four 
dimensional Lichnerowicz operator plus a d 2 piece. Performing a fourier decomposition of 
h ab m the fifth dimension yields 

A L h ab = (A 4 - m 2 )h ab = 0. (1.2) 

Since the four dimensional Schwarzschild Lichnerowicz operator has no unstable modes, 
adding a mass should only increase stability, hence it has been conjectured black strings 
are stable. 

On the other hand, glibly speaking, horizons are like soap bubbles, they have a surface 
tension, k - the surface gravity, and soap bubbles do not like being cylindrical! More 
formally, a portion of horizon of length L contains mass M. = ML and has an entropy 
proportional to M 2 /L. A five dimensional black hole on the other hand has entropy 
proportional to M. 3 / 2 . Thus for large lengths of horizon, the mass contained within the 
horizon contributes a much lower entropy than if it were in a hyperspherical black hole. 
This indicates that for large wavelength perturbations in the fifth dimension, we might 
expect an instability. 
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As it turns out, the thermodynamic argument is correct. The flaw in the other ar- 
gument is to assume that the number of physical degrees of freedom of the effective four 
dimensional tensor gauge field, h a b, remains the same as in Einstein theory. Clearly this is 
not so, for the fifh dimension adds an effective mass to h a b. A massive tensor field in four 
dimensions has three degrees of freedom as compared with only two for a massless one. 
We will show that this extra degree of freedom has a spherically symmetric mode which is 
responsible for the instability. 

The stability of the five dimensional black string was investigated analytically 9 , with 
the conclusion that there is no non-singular single unstable mode on a Schwarzchild 
time t = surface, however, this argument did not prove stability. As emphasized by 
Vishveshwara 10 in his original Schwarzschild stability argument, the non-existence of a 
single unstable mode does not preclude the existence of a composite unstable mode, with 
the combination cancelling the singular behaviour of an inadmissable single singular mode. 
An example of this sort of situation occurs in the coloured black hole instability, recently 
confirmed by Bizon and Wald 11 . That this was probably the case for five dimensional black 
strings was first pointed out by Whitt 12 , who analyzed the stability of Schwarzschild in 
four dimensional fourth order gravity - a different physical situation, but mathematically 
identical equations to those studied in ref [9]. The key difference in Whitt's analysis was 
the choice of an initial data surface ending on the future horizon (see Figure 1.). By avoid- 
ing the neck of the Schwarzschild wormhole, one avoids the fixed point of the isometries 
used to generate the mode decomposition, which avoids in this case issues of superposi- 
tion. In any case, such a choice of initial data surface is mandated in any real physical 
scenario, since a black hole or brane must form from collapse, and hence will not possess 
a Schwarzschild wormhole. 

In a previous letter 13 , we showed that using an appropriate initial data surface, un- 
charged black p-branes were unstable, and conjectured that the same would hold true for 
charged black p-branes. In this paper, we extend the analysis of [13] to include charge, 
and show that the instability does indeed persist. It is worth stressing that this instability 
is not of the Reissner-Nordstrom form - hidden behind the event horizon, but it is a real 
physical instability of the exterior spacetime which could potentially fragment the horizon. 
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It is important to emphasize that this can occur classically, for although under regular 
conditions horizons do not bifurcate 14 , if one has a naked singularity, then bifurcation 
is possible. Since an instability calculation is by its nature linear, it cannot predict the 
endpoint of an unstable evolution. However, the entropy argument does lend support to 
the fragmentation scenario and violation of the Cosmic Censorship Hypothesis. 

The layout of the paper is as follows. In section two we summarise the black p-branes 
we are investigating, namely those of Horowitz and Strominger 7 , and derive the general 
perturbation equations. We discuss gauge constraints and the issue of boundary conditions 
in detail, then we set up the mode decomposition for our analysis. In section three we 
specialise to zero charge and show how this greatly simplifies the problem, and present an 
instability. We then use this 'zeroth order' solution to generate an instability for small 
charge in section four. In section five we present the general case, for arbitrary charge, 
and present our conclusions in section six. 



2) The perturbation problem. 

We start this section on the derivation of the perturbation equations by reviewing the 
particular black p-branes we will be investigating. We focus on solutions of low energy 
string gravity with symmetry IR 10_i:, xIR t xS D ~ 2 and a possible 'magnetic' charge. The 
action for such a theory is 

J d w xV^9e-^[R + 4(V0) 2 - j^yF 2 ] (2.1) 

where F is a (D — 2)-form field strength. Varying this action yields the equations of motion 

□^-2(V^) 2 + F 2 ^f| [ = (2.2a) 

V ai [e- 2 ^F ai -- aD ~ 2 ] = (2.26) 
R ab + 2V a V b </> - {D 2 _ 3)l F aa2 ... aD _ 2 F b a2 ---- aD - 2 = (2.2c) 

Searching for a solution with the required symmetry indicates a metric of the form 

ds 2 = -e A dt 2 + e F dr 2 + e B dx l dx t + C 2 dQ 2 D _ 2 (2.3) 
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where all metric functions are functions of the radial variable only, and the index i runs 
from D + l to 10, and represents the coordinates in the p-plane of symmetry of the p-brane. 
Looking for a 'magnetically' charged brane, 



F = Qe D - 2 (2.4) 

where is the area form on a unit (D— 2)-sphere, it is straightforward to show that a 

solution exists and takes the form 7 



K r 1 ' (2.5) 

C(r)=r , e-^ = l-C-) D -* , Q 2 = ^(r + r_)^" 3 

If r_ = 0, then Q 2 = = 0, and the solution is the uncharged Sch^ xIR 10-15 , which we 
have argued to be unstable. However, in the presence of charge there is also a non-trivial 
dilaton field, and it is unclear how this will affect the physical instability. We first vary 
the equations of motion to obtain the bare perturbation equations, then discuss bound- 
ary conditions on our perturbations in terms of a non-singular (or generalised Kruskal) 
coordinate system on the event horizon, then return to the perturbation equations, dis- 
cussing gauge choices to simplify these, finally discussing an additional simplification to 
the electromagnetic perturbation. 

In order to write the perturbation equations, we use the usual notation 

Sg a b = h ab (2.6) 

and 

hab = h ab - \hg ah (2.7) 
varying the equations of motion can be seen to give 

V ai [e- 2(p (5F) ai - aD - 2 ] - 2(V ai 5(p)e- 2(p F ai - aD - 2 — (D — 3) e - 2 *F aib[ - aD - 2 V ai hl 2] 

-h aib V ai [e- 2 ^F b a2 - aD - 2 } - V ai (^)e" 2 ^- fflD - 2 =0 (2.8a) 
nScf) - h ab V a V b (p + 2h ab V a (t)V h (t) - V c (pV b h bc - 4V a <W> 
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- ^^F ca2 ... aD _ 2 F d a ^ a ^ + 2^-^dF ai ... aD _ 2 F^--^ = (2.86) 
Dh ab + 2R cadb h cd - 2R e{a h e b) - 2V (a V\ e \h e b) - AVaVbdcj) - 2V c (j)V c h ab + AV c (jN {h h c a) 
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+ 



2F {a a *- a °->6F b)a2 ... aD _ 2 — (D — 3)h cd F aca3 ... aD _ 2 F bd a ^ a ^] = (2.8c) 



(D-3)\ 

These are the bare perturbation equations, we must now examine boundary conditions, 
and whether there is any choice of gauge which simplifies (2.8). 
Boundary conditions. 

Now we turn to the question of boundary conditions, which are a key to this problem. 
Obviously, we want to place initial data on a Cauchy surface for the exterior spacetime, but 
such a surface necessarily touches the horizon, which is singular in Schwarzschild coordi- 
nates. There are therefore two issues here: One is how to define 'small' for the perturbation 
at the horizon, and secondly, which initial data surface to impose these constraints upon. 

The first issue is straightforwardly dealt with. Although the horizon is singular in 
Schwarzschild coordinates, it is not a physical singularity, merely a coordinate singularity. 
In four dimensions, non-singular coordinates have been known for some time - Kruskal 
coordinates. In order to generalise these to higher dimensions, the starting point is to 
identify the generalised tortoise coordinate 

dr 



e (F-A)/2 dr = / r = _J_ 



1-Kin 



n =i V ~ 6J (l-e°-*r+/r) 
— — — r+log(r-e^-3r+) (2.9) 

n=l ^ ' 

So, for example, 

r\ = r + r + log(r - r+) 

r * = r + \r+ log (2 .10) 
r* = r + \r + log(r - r+) - ^ log(r 2 + rr + + r\ ) + tan" 1 [^^j 
etc. In particular, as r — > r + 

r* D -^^\og(r-r + ) (2.11) 
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Having the tortoise coordinate, we can now define the Kruskal null coordinates 



P + =exp{(D-3)(t + r*)/2r+} 
P_ = -exp{-(D-3)(t-r*)/2r + } 

In terms of which the radial-time part of the metric looks like 



(2.12) 



, 4r 2 (D-3) 2 e A 
ds\ = + p+p _ dP+dP_ (2.13) 



which is finite as r — > r + . 
Finally, we set 



R = P+-P. 
T = P + + P_ 



(2.14) 



in terms of which 



ds l = - r +( D ^^ (dR 2 - dT 2 ) = f 2 (r)(dR 2 - dT 2 ) (2.15) 
P+P- 

in order to readily identify our initial data surface. 

This leaves us with the problem of an initial data surface. The domain of dependence 
must obviously include X + , thus a surface touching the future horizon, or the neck of the 
Schwarzschild wormhole is acceptable, but a surface touching the past horizon is not, unless 
it passes through and extends to the opposite horizon on the Penrose diagram. Following 
Whitt, we impose regularity in the Kruskal system as r — > r + at constant non-zero T, i.e.,, 
R — > T from above as depicted on Figure 1. This avoids the issue of mode superposition 
discussed earlier, and secondly, we believe it to be a better physically motivated choice 
of surface. This is because in practice a black hole (or brane) would form in a collapse 
situation, and hence would not have a Schwarzschild wormhole; analyzing stability would 
necessarily require a surface ending on a future event horizon. 

In terms of the Schwarzschild components 



TT (D - 3) 
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h 11 = 



f 2 (r)(T 2 -R 2 ) 



R 2 h u + 2RTh^ + T 2 h^ 



(2.16a) 



h RR = 



h TR = 



h T * = 



h R * = 



[D-sy 



f 2 (r)(T 2 -R 2 ) 

(D - 3) 4 
~P{r){T*-R*) 

(-P-3) 2 
f(r)y/(R*-T*) 

(-P-3) 2 
/(r)#^j 



T 2 h U + 2RTh U + R 2 h^ 
RT(h U + h**) + (T 2 + R 2 )h^ 
Rh b + Tff* 
Th { * + Rh™ I 



(2.166) 
(2.16c) 
(2.16d) 
(2.16e) 



all of which should be finite as r — > r + on our initial data surface. Since = T + 0(r — r + ) 
as R — > T from above, this implies that the various combinations of the normalised /i a6 
within the square brackets above must be 0(r — r+) as we approach the event horizon. 
Now we turn to simplifying the perturbation equations by a judicious choice of gauge. 
Gauge considerations. 

In gravity, physics must be invariant under a re-labelling of coordinates. Such gen- 
eral coordinate transformations (get's), are generated by vector fields £ a , the effect of an 
infinitesimal get being to push the coordinates an infinitesimal amount along the inte- 
gral curves of £ a , such that x a — > x a + £ a . Under such a gauge transformation, physical 
quantities transform as 

V -> L(P (2.17) 

where denotes the Lie derivative with respect to £. Therefore a pure gauge perturbation 
is of the form 

h^ab = 2 £( ;6) 

8<f> = i a V a <t> (2.18) 

&Fa\...a D -2 = F ai ...a D - 2 ;d£ + (D — < ^)Fd[a 2 ...an-2^;a 1 ] 

It is straightforward to verify that such a pure gauge perturbation satisfies our generalised 
Lichnerowicz perturbation equations. 

Using such a pure gauge perturbation, we see that 



V a hf = R b J a + n£ 6 
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(2.19) 



However, the right hand side of this equation is the curved space wave operator, which 
given a source and initial conditions can be integrated to give a solution for £ a throughout 
the manifold 15 . Thus, if the divergence of h a i is not initially zero, we can change gauge in 
order to make it so. This simplifies (2.8) to the following: 

V ai [e- 2 ^(5F) ai - aD - 2 ] - 2(V ai 5(/))e- 2 *F ai - aD - 2 - (D - 3)e~ 2(l, F ai ^- aD - 2 V ai h a b 2] 

- h aib V ai [ e -2<t> Fb a2 - a °- 2 ] = o (2.20a) 
DS<j) - h ab V a V b <p + 2h ab V a (tN b (t) - 4V a oc/>V> 

~ J^^F ca2 ... aD _ 2 F d a ^ a ^ + 2^f| [ oF ai ... aD _ 2 F--— = (2.206) 
Uh ah + 2R cadb h cd - 2R e(a h e b) - AV a V h 5<j) - 2V c cj)V c h ah + AV c cjN {b h c a) 



+ 



2F (a a *- a °- 2 5F b)a2 ... aD _ 2 — (D — 3)h cd F aca3 ... aD _ 2 F bd a *- aD -> = (2.20c) 



(D-sy. 

The residual gauge freedom are those get's which satisfy 

R b a C + ae = o (2.21) 

Notice that in a vacuum spacetime, i.e. no dilaton, no 'electromagnetism', the trace of the 
metric perturbation equation becomes 

nh = (2.22) 

and we can choose an harmonic £ a to set h = 0. Therefore in vacuo we may make the 
additional gauge choice of trace-free. However, more generally it should be noted that 
tranversality is all that can be assumed. 

Since a pure gauge perturbation automatically satisfies the perturbation equations, 
deciding whether a putative instability is physical or not reduces to investigating whether 
it can be expressed in this form. Alternatively, one can identify physical perturations by 
using gauge invariant variables. It is straightforward to adopt Bardeen's results 16 to the 
case at hand. The idea is to choose a particular linear combination of variables and their 
derivatives such that the result becomes independent of £ a . The details and a particular 
set of gauge independent variables are given in Appendix A. 
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Finally, before beginning the analysis of the equations in detail we will make a further 
simplification by showing that the 'electromagnetic' perturbation, 5F a b, can be taken to be 
zero. In order to perform a stability analysis we will Fourier decompose all the perturba- 
tions in terms of the symmetries of the background spacetime. Since Kruskal coordinates 
are not well suited to such a decomposition, we perform our analysis in Schwarzschild co- 
ordinates, transforming to Kruskal at the horizon to check boundary conditions according 
to (2.16). The Fourier modes in the time and p-brane directions are easy to identify, for an 
instability they are of the form e,^ t+lfJ,iX S ph er j ca i harmonic modes will depend on 

the number of dimensions, D, that the black hole sits in, as well as on the tensorial nature 
(scalar, vector, and so on) of the perturbation being analysed. However, since in general 
higher angular momentum modes are more stable, and since our advertised instability is 
an s-mode, zero angular momentum modes are all we shall be considering. This clearly 
means we can make no stability claims if we fail to find an instability, however, since our 
primary concern is to chart regions of instability, an s-wave analysis will suffice. 

Therefore, the spherically symmetric perturbations take the form: 



S(j) = e nt+iliixi f(r) 



5F 



q{r)cD-2 angular components 

mixed angular / t-r-i components 

8Ft t-r-i components 



(2.23a) 
(2.236) 



-iftf(r) 












H tt R tr 








H rj 
















K(r) 














Kj sin 2 6 


only if D 


-2 < 5. 







(2.23c) 



We may now see that the F-perturbation equation, (2.20b), becomes simplified, for 
noting that the background field F = Qeo-2 has only angular components, and that 
hp = K(r)5p/r 2 for angular components, 



F aib ^- aD - 2 V ai hl 2] = F ai ^- aD - 2 h a h 2 l +F aib t- aD ~ 2 T a r 2 n \h c h = 



b,ai 



cai lb b 
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and 

V r/t [e- 2 ^F a ; t - aD - 2 ] = h a ^ b V ai [e-^F^- aD - 2 ] oc g ab V ai [ e -^ F «-— ] = 

Hence the F-perturbation equation reduces to 

V ai [e- 2 HSF) a b 1 - aD ~ 2 ]=0 (2.24) 

We can make the further observation that since F and 5F are derived from a potential, 
d(6F) = 0, hence 

5F 0<j) ,t = Slq{r)ee<t, = (2.25) 

provided Vt ^ Q. Thus only 5Ft can be non-zero, and this is independent of the angular 
variables. 

We may now argue that 5F T also vanishes. Without loss of generality, we will show 
this for D = 4, the generalisation to higher dimensions being an iterative process. 
Consider d(SF) = 0. 

F a b,c + Fbc, a + F caib = (2.26) 

For a = t, b = r, c = i, we have 

8F' Tit + n5F T „ + imSF Ttr = (2.27) 
Setting f t = T<(ii5F T i t and f r = T,fXiSF Tir and summing (2.27) appropriately gives 

f' t - Qf r + ifi 2 5F Ttr = (2.28) 
where fx 2 = T,fi 2 . Now, we return to the equation of motion (2.24) which implies 

QSF*? + i^F l T a + 5F£ a ' + 2 -bF r T a = (2.29) 



For a=r, we have 

(r — r_) 



(r-r+) 

12 



n5F Ttr + if r = (2.30) 



Substituting from (2.28) implies 



fl = 



\j?{r — r+) 
0(r — r_) 



+ 



fr 



(2.31) 



Whereas taking (2.29) for a = i, multiplying by [ii and summing gives 



-tog" ft - (g rr f r )' - -g rr f r = 
r 



Substituting for f r from (2.31), and rearranging gives finally 



(2.32) 



(r-r-)(r-r + ) ( (2r - r+ - r_) 
I2 ~+~ 



-l„ 2 + r> 2(r " r ~ 



(r — r+) 



/ u 2 (r_|_ — r_)(r — r + ) 
r[// 2 (r — r_|_) + 2 (r — r_ 

/t = 



(2.33) 



By inspection of this equation we can see that there are no regular solutions and thus we 
must take f t = 0. Hence f r = 5F rtr = 0. 
Now, returning to (2.27) implies 



SFUt + ^Kn = 



(2.34) 



and (2.26) for a = i, b = t, c = j implies 



iHjSF Tit + i(jndF T tj + Q5F Tjl = 



(2.35) 



Multiplying by fii and summing gives 



i/i 2 SF Ttj - Qfj = 



(2.36) 



where fj = T,HidF Ti j. Taking the j-component of the equation of motion 



nSF Ttj g u + if j + (SF Trj g rr Y + -SF Trj g rr = 0, 



substituting in for fj , SF Tr j, and rearranging gives 



(r ~r-)(r -r+) (2r-r + -r_) 
12 bh T tj + -3 d *Ttj 



(2.37) 



-^ + n 2 ^—^SF Ttj =0 (2.38) 
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As before, this has no regular solutions, hence SF rt j = 0. (2.34) and (2.35) then imply 
that all other components of 5 Ft vanish, hence the perturbation 5F ab may be taken to be 
zero in all future analysis of the spherically symmetric perturbation equations. 

Finally, to summarize, the equations of motion for the spherically symmetric pertur- 
bations we are interested in analysing reduce to the following 

nScf) - h ab V a V b (f) + 2h ab V a <pV b <p - AVJ(t)V a (t) 

h cd 



a,2---a.D — 2 



^jj _ ca 2 ...a c _2- t d 



= (2.39a) 



Dh ab + 2R cadb h cd - 2R e{a h e b) - 4V a V 6 50 - 2V c (fN c h ab + AV c <jN {b h c a) 
~ J U ^h cd F aca3 ... aD _ 2 F b ; 3 --- a ^ = (2.3%) 

the 'electromagnetic' perturbation vanishing in this case. Replacing the background values 
of the dilaton and Riemann tensor lead to the final versions given in appendix B. 

3) The zero charge case. 

We will first give a detailed analysis for the case when the black string or brane is 
uncharged. This simplifies the calculation for two reasons. The first is that the form of the 
metric is much simpler, leading to more tractable coefficients in the perturbation equations 
of Appendix B. The second is that as can be seen from appendix B, some equations decouple 
from each other; in particular, the perturbation of the dilaton can be neglected and we 
can deal solely with the metric perturbations. In addition, as we have mentioned before, 
it is consistent to take a vanishing trace of the metric perturbation. 

We can first investigate the transverse terms of the metric perturbations in the 10 — D 
dimensions. These perturbations transform like scalars under the D— dimensional coor- 
dinate transformations as is well known in Kaluza-Klein decompositions. They are com- 
pletely decoupled from all other perturbations and obey equations of the form 

Ah u " + Bh"' + Ch" = (3.1) 
where A, B and C are function of r. Their explicit form can be found in Appendix B. We 
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can investigate the behavior of h n as r — > oo and and as r — > r+. We obtain 



/i" ~ exp + /U 2 r — > oo 

(3.2) 

~ (r - r+) ±nr + r -> r+ 

The regular solutions correspond to the minus root at infinity and the positive root at 
the horizon. It is therefore easy to see that a necessary condition to obtain a regular 
solution is to have a change in sign in the ratio of coefficient A/C in equation (3.1). From 
appendix B it is seen that this ratio does not change sign for the transverse part of the 
metric perturbations. We therefore set these perturbations to zero in our search for an 
instability. 

We can now study the vector perturbations W 1 . An analysis similar to the scalar case 
shows that they do not lead to an instability, they are therefore set to zero. We are thus 
left with the tensor perturbations. 

The tensor perturbations are more complex as the equations are coupled with each 
other. It is possible to rewrite them in term of a single variable, h tr say, using the gauge 
conditions described in appendix B: 

o = - + (D - 3) [j? 



rV 



, Q 2 [(D-2) + (2D-7)( r -±) D ' 3 ] _ 3(D-3f( r -±) 2{D - 3) [(D-2)-( r -±) D - 3 ] ^ Htr/ 



rV Ar 3 V 



+ { ^ + nyv) 2 + " 2[4( "- 2) - 8( "- 2) ^)7;-,f- 34 " +5 " 2) ^) 2(g " 3)] 

^[ 4 (D-2)-4(3D-7)(^) g - 3 + (^+2D-ll)(^) 2( °- 3) ] 

( D -3) 2 (^)^- 3 )[(D-2)(2D-5)-(D-l)(D-2)(I±) g - 3 + (I±) 2(g - 3) ] 1 
"I" 4r 4 V 2 \n 



(3.3) 



Before dwelling on the existence of an instability for the tensor perturbations we would 
like to comment on the flat space case, i.e r + = 0. In this case the ratio of the coefficients 
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of the first and last term does not change sign and we get the expected result that there 
is no instability. 

When O > (D — 3)/r+ it is possible to use the same analytic argument as for the 
transverse perturbations to show that no instability exists. However when O < (D — 3)/r+ 
the ratio of the coefficients of the first and last term does indeed change sign and we must 
resort to numerical techniques to establish the existence or otherwise of an instability. 

It is possible to reduce the set of second order equations to first order ones by using 
the two gauge conditions (B.lOa-b). By differentiating these and using the second order 
equation we reduce them to first order equations but as we have only three variables 
and four first order equations we can reduce one of them to a constraint between the 
fields themselves (without derivatives). One of the variables can thus be reexpressed as a 
function of the others. 

First let us define 

zrrr 

H± =V+H U ± — - (3.4a) 
H = -H tr (3.46) 

where V+ = (1 — (r + /r) D ~ 3 ). The Lichnerowicz and the gauge equations reduces to the 
constraint: 



H-[^ + ^± + { -^(-(D-3) + (D-4)V + + Vf)] 
+ H?±[-Cl(D-2)-£{(D-3)-(D-l)V + )] 



(3.5) 



and the following two first order differential equations 



and 



^ = JL( g+ + g _)-K J -»)+^] g (3.6) 
dr 2VV + ; rV+ V ; 



dr fl 2r + 2rV+ V ; 
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^From this we can deduce that as r — > oo 



H = ^Q? + ^F+e r ^ n2+ » 2 - x/O 2 + //2 j p_ e -r-V n2 +/^ 2 (3.8a) 



and as r — > r+ 



if = (D - 3)rf " 4 (-i + ^)G + (^" 3 - rf" 3 )"^ 

— (D — 3)rf " 4 (-i - ^)G_(r^ 3 - r^" 3 )" 1 ^ (3.9a) 
if_=(^ + ^)G + (r^- 3 -rf- 3 )S 

S 2 r_|_ 

The behavior of can be obtained from the constraint (3.5). 

The regular solution corresponds to F + = and G- = 0. We have used a Runge-Kutta 
algorithm with variable stepsize to integrate these equation from r\ = 200 to r 2 = 2.000002 
with a tolerance e = 10~ 6 . We do this integration for fixed ji and various values of O. Near 
r = r+ we calculate the ratio R = G-/G+ and look for a change in sign (as O varies). 
When a change in sign occurs we home in towards the value of O for which this change 
occurs and ensure that the ratio decreases towards such a value (an increase would imply 
that G- goes through a zero). Figure 2 shows the behavior of the function H u , H tr and 
H rr . In Figure 3. we have plotted Q as a function of \i for which an unstable mode has 
been found. This plot takes r + = 2, but behavior for other values of r + can be obtained 
from the scale transformation 

r+ — > ar+ 



~ a (3-10) 



There are at least two interesting points in Fig. 3. The first one is to note that O does 
not seem to go to zero as \x — > 0. This might be a numerical artifact as our boundary 
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conditions are different in this limit. However if it is not a numerical artifact it would be 
surprising at first, as in this particular case the equations reduce exactly to the stability 
of D-dimensional black holes. We know that black holes in four dimensions are stable, 
moreover the mode under investigation here is the spherically symmetric mode and we 
know by Birkhoff 's theorem (and its generalization to higher dimensions) that the solution 
must be Schwarzschild. We must therefore conclude that this mode is pure gauge. We can 
in fact work out what this gauge transformation looks like in the asymptotic regions. 

As the system remains spherically symmetric and has no dependence on the extra 
dimension, we must have 

£° = (£*,f\0,..,0) (3.11) 
It is easy to verify that the particular gauge (in the large r limit) 

= f = e ^ f - r) (3.12) 

corresponding to metric perturbations of the form: 

h tt = £f,t „ _Q e «(t-r) (3.13a) 
h rr = ^r;r ^ _ ne n(t-r) (3.13b) 

h tr = £(*;»•) « -Q e n(t - r) . (3.13c) 

This gives rise to a transverse and traceless metric as can be easily verified. In the r — > r + 
limit we have 

e = -r+(r - r+)- 1+nr + ; f = (r - r+)^+ (3.14) 
corresponding to the metric perturbation 

h u « {-\ + Or+)r+(r - r+ )- 2 +^+ e n * (3.15a) 

h rr w ( ~ 1 + 2 " r+) ( r - r+) nr +e nt (3.156) 
h tr w -(-- + Or+)(r - r+)- 1+nr +e n * (3.15c) 

Of course we still have to show that there is a regular solution which links these two 
asymptotic behavior. This can be done numerically. 
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The result for \i = might leave the reader worried that the instability found might 
be pure gauge. However the case where \x 7^ is different. In this case it is easy to convince 
oneself that they cannot be pure gauge. The simplest way is to assume we can write these 
perturbations in term of £ a and find a contradiction. Because perturbations in the 10 — D 
dimensions are zero, we would have 

h<A J o j i) • i io 

thus this implies that £ J = for /jj ^ as h?i = 0. If we now look at = this 
equation implies that £ M = and thus h^ v = 0, a contradiction. This ensures that the 
perturbations found in the uncharged case cannot be pure gauge and therefore must be 
physical. In the charged case the perturbations in the extra dimensions are non-zero and 
thus such a proof is more difficult. In this latter case we will make use of gauge invariant 
variables described in appendix A to show that the perturbations are not gauge artifacts. 

It is interesting to understand where this extra degree of freedom for the gravitational 
field has come from. It is a result of the degree of freedom of the gravitons due to the 
extra dimensions. One easy way to make this clear is to use a Kaluza-Klein analogy. 
The gauge transformations in D dimensions give the graviton degrees of freedom which 
are transverse to the direction of propagation. This is related to the masslessness of the 
graviton. Once the extra dimensions come into play, through the dependence exp i[ijZ^ , this 
effectively gives a mass to the graviton and thus a longitudinal component corresponding to 
a spherically symmetric mode in D dimensions. It is this component which is responsible 
for the instability. In the limit as fx — > 0, this mode becomes pure gauge as the mass 
disappear. 

The other interesting point to notice in Fig. 3 is the fact that there is a maximum 
value [i™- ax for the instability to exist. The analysis carried above can be repeated in the 
case where the extra dimensions are periodically identified. The only difference in this 
case is that the values of \ij will be quantized as 2nn/Lj, for n = 0,1, ... and L being the 
length of extra dimensions. ^From this we can see that if the Lj are small enough, the first 
value of Hj will be larger than \x™ ax and thus the instability will not occur. The instability 
occurs only when the black holes (in D dimensions) have a Schwarzschild radius of the 
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order of the size of the extra dimensions as suggested by the entropy argument described 
in the introduction. Thus if we believe the assumptions of string theory that the 10 — D 
extra dimensions are wrapped in small circles, the instability found here will not have any 
effect on astrophysical black holes. It is interesting that this instability appears to break 
the duality of the solutions obtained by replacing the radius L of the extra dimension by 
its inverse 1/L. If we choose L appropriately, there is an instability for this spacetime 
but there will not be one when L — > 1/L, however, as we have considered only momentum 
modes in this calculation, and not winding modes, it is possible that there will be a winding 
mode instability for the dual spacetime. However, we have neither verified nor refuted this 
claim. 

We have discussed the existence of the instability but what does it correspond to? All 
the perturbations have the form expi^jz^ and thus implies an oscillatory behavior of the 
distance scale as a function of the extra dimensions. To make this clearer we can work out 
the location of the apparent horizon. To do so we introduce ingoing Eddington-Finkelstein 
coordinates, and for ease of description restrict our attention to D = 4: 

ds 2 = - ( r ~ r+ ' du 2 + Idudr + r 2 dfl 2 + dz 2 (3.16) 
r 

where u = t + r + r + ln(r — r + ) in terms of the usual Schwarzschild coordinates. The 
apparent horizon is located where the outgoing light-cones have zero divergence. The 
spherically symmetric light rays for this metric are given by 

u = const. ingoing 

dr _ (r-r+) .„ • „ (3.17) 



AT = —2T 1 outgoing 

In the case of the unperturbed black strings the apparent horizon is at r = r + . If we 
perturb the metric slightly we get for the new outgoing null geodesies 

huu ~ — r +2 d^ {1 + hur)+hrr \d^) = ° 
which to first order gives 

dr_ (r-r+) (r - r+) 1 (r - r+) 2 
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which has a zero for 



r & r + (l + h uu ) + lim (r - r + )h ur + 




h 



rr 



(3.19) 



Using the behavior of the solution obtained above we finally get 



r ~ r + + const. cos(/iz) 



(3.20) 



This shows that the apparent horizon of the perturbed spacetime oscillates as a function of 
the extra dimensions. The schematic behavior of the apparent horizon is depicted in Figure 
4. Thus it would appear that this perturbation destabilises the event horizon causing it to 
ripple in the transverse dimensions. 

4) The small charge case. 

We have seen in the previous section that uncharged black strings and p-branes are 
unstable. We now investigate the question whether a charge on these black objects might 
stabilize them. Physically, it seems reasonable that they will also be unstable at least for 
small charges, this is because for small charge the metric is essentially the same as the 
uncharged case. As long as we remain outside the event horizon (r > r+), as indeed we 
do to investigate the instability, the effect of the charge should be negligible. However 
a non-trivial effect of the charge is that now the 10 — D dimensional perturbations that 
previously vanished become coupled to the D-dimensional non-zero perturbations. This 
now complicates the issue as we can no longer set the former perturbations to zero; all 
the perturbations must be solved simultaneously. An additional complication is that it is 
not possible to take the trace h = g a bh ab to be zero, for as we have already shown, this 
simplification is only possible in vacuo. The perturbation of the scalar field 4> couples to h 
and prevents us from setting that part of the metric to zero. 

In order to guide the reader through the maze of these equations we will first focus on 
the small charge case. Indeed in this case the problem, although more complicated than the 
uncharged case, is drastically simplified from the general charged case. The simplification 
come through the observation that the coupling of D-dimensional perturbations to the 
10 — D ones are always through a factor proportional to the charge Q 2 = r_r + /2. As the 
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10 — D dimensional perturbations are vanishing for Q = 0, we deduce that they are in fact 
proportional to Q at small charge. This lead us to an expansion of the perturbations in 
powers of r_. Using such an expansion will be zeroth order, and the other variables, 
h?^ , h? k and 6(f) will be first order. Notice that h, the trace of h ah , will also be first order. 

Before starting the analysis, note that from the form of the perturbation equations in 
Appendix B, it is easy to see that writing 

hP* = —h» z (4.16) 

removes the (10 — D) z-degrees of freedom to just one transverse degree of freedom, which 
we call l z\ Additionally, for calculational simplicity, we will consider the specific case of 
the black 6-brane (D = 4), which will at least remove the variable D from our equations! 
The analysis can be extended to other black strings and p-branes, however, since we have 
already shown black p-branes to be unstable for all values of D in the previous section, an 
analysis for general D at this stage would be neither economical nor illuminating. 

We first analyze the equation for the perturbations of scalar field /, h zz and the trace 

h: 

+ - U + J^r) J + r - <r ' 2 2r+> A- 

r r z \ (r — r + ) J r A 

r-(4r-5r + ) ^ rr _ r_r + (r-r + ) feft = Q 
4r 3 (— r + r + ) 4r 5 



(4.2) 



(r T+) h zz " + (2r o r+) h zz ' -(v 2 + , ^ h zz + 2fi 2 f(r) = (4.3) 
r r z \ (r — r+) J 

^^-V")+ { ^<.Uf)-U + j^) (l-2f)- 2 -^±K = (4.4) 
r 2 r z 2 \ {r — r + )J 2 r z 

There is also an equation for h rz and h tz and the gauge condition 

V a h az w h rz ' + i/jh zz + nh tz + -h rz -i^h^O (4.5) 

r 2 
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Taking a derivative of this gauge condition and using the h tz equation we get 



G^W' - if ft' + »'■'+#/') + (f + ^ry) ft" + ff ft" + |±ft" - o 

(4.6) 

This gives us two first order equations for h rz and h tz 

With these equations we can investigate the behavior of the fields at infinity and near 
r = r + . In the large r limit all these fields have a behavior of the form e ±r v /n2 + At2 , the 
regular solution having the negative sign. The coefficients of /, h, h zz , h tz are independent 
of each other but once these are fixed, the coefficient in front of the exponential of h\ z 
(denoted by a subscript i) is given by 



hV 



^==^ hi ^-, hr ). (4.7) 



in terms of the coefficients of the exponential of the other perturbations. Near r = r + , the 
perturbations h, /, h zz , h rz behave as (r — r+) ±fir+ and h tz as (r — r + ) _1±slr +. 

With these boundary conditions we can integrate eq(4.2-6) to find the unstable mode 
assuming the zeroth order solutions found previously. We have used a similar technique 
as in the uncharged case to find the regular solution near r = r+. By adjusting the values 
at large r we have been able to find regular solutions for all the perturbations. Fig. 5 gives 
the behavior of the modes which lead to the instability. 

5) The charged case. 

We have now shown that black strings and branes with small charges are unstable 
by assuming a perturbative analysis in term of the charge. As we have seen in that case, 
we have a zeroth order equation in a small charge perturbation for the four dimensional 
metric perturbation. The dilaton and the extra dimensions are first order and could be 
solved assuming the zeroth order solutions for the four dimensional metric. In this section 
we give the argument for the general charge. The idea of how to handle the equations is 
very much the same as in the previous cases but the calculations are much more tedious 
as the equations will now involve the charge parameter (or rather r_ = 2Q 2 /r+). 
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The equations of motion for the different fields are given in appendix B. As before, 
we were unable to show the existence or absence of an unstable mode analytically. We 
have thus resorted to numerical evaluation of the unstable mode. We can rewrite the 
set of equations in Appendix B as a set of first order equations suitable for numerical 
integration. First let us define a new variable q = h/2 — 2/ where h is the trace of the 
metric perturbation and rewrite the equations for /, q, h zz as 



7T f = - 



(r — r_)(r — r + )\ (r — r + ) 

r_(4r 2 — 5rr_ — 5r + + 6r_r + ) 



2 n 2(r-r-) r-(r-r + )(r--r + ) /tft 



4r 3 (r — 



r_(r — 2r+) (2r 2 — rr + — r_r + ) 



4r 2 (r — r_) 2 (r — r+) 
-K + 



TTf 



7V 



i^ ± K + { l ?+^ 7 -—^)q 



(r — r_) 



q (r — r_)(r — r+) V r 2 
(r — r_)(2r — r + ) 



(r-r+) 



1.2:2' _ 



n hzz — 



(r — r_)(r — r+) 



(- V/ + 



i2\xr- 
r(r — r_) 



{ji 2 + s 



(r — r_) (— 2r + r_+r + ) 



(r — r+) 



)/i Z2 + 



(5.1) 



The next step is to use the gauge condition (B.lOd), the gauge condition for the index 
L z\ The equation for h tz ' is obtained by taking the derivative of this and using the equation 
of motion of h rz to eliminate the second derivative. This gives an equation in terms of the 
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variables and their first derivatives. Explicitly we have 



If 0(r_-r + ) tz firr. 



ih lz = - — v ~ -LL ih lz 1 h 

fi\ (r — r-)(r — r+) (r — r_) 2 (r — r+) 

+ ( (r _ r+) 2 + ( r _ r _)( r _ r+ ) )^ r2 - 2 W + ^ " ^ hZZ ) ( 5 - 2 ) 
ifr™' = - m - 2/lf - Qih tZ - ( 2r ~ 3r ~) h rz + h »» 

ryr — r_) 

We can use the gauge condition for the indices t and r to get first order equations for h tr 
and h rr 

h ir< = _ n(r-r.)( q + 2f) _ nhtt _ (2r*-4rr--rr + + 3r r + ) htr _ u 
2(r — r + ) r(r — r_)(r — r + ) 



(5.3) 



hrr , _ (r - r-){r - r+) + . _ (r - r + )(-2r 2 + rr_ + 3rr + - 2r_r + ) feff 
^.3 • 2r^(r — r_) 

- - (6r2 - 9 ; r -- 7 :; ++1 ° r - r+) ^ - 

2r(r — r_)(r — r+J 

^ " + ~2 K + 27r /) 

Taking a derivative of the gauge condition with index t and getting rid of the second order 
derivative using the second order equation for h tr we get 

At' _ 1 ( °( r - - r +) „ °( r - - r +) z,tt 



/i" = i_ —--q — — — h + 

fiv (r — r + ) 2 (r — r-)(r — r+) 

/ ° 2r2 + A 2 x fofr + ^(r- - r+) | 

(r — r + ) 2 (r — r_)(r — r+) (r — r_)(r — r+) 

, fi 2 rr- ^ Or(-2rr_ + rr + + r_r + ). rr 

0(r — r_) 2 (r — r + ) (r — r_)(r — r_|_) 3 

- ^(r^ + 7 ^ - (£ - - 2-/) 

O v (r-r + ) 2 (r - r_)(r - r + ) y v fi (r - r+) 7 



r 2 

+ 
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We have now the field equations written in terms of first derivatives but we still have 
equation left which comes from the time derivative of the gauge condition with index 
This leaves us with a constraint between the variables in equations (5.1-4) given by 



20 2 (r r 2 2jU 2 (r — r_)(r — r + ) (r_ — r + )(— 2r 2 + rr_ + 3rr + — 2r_r + ) 



+ ( 



8^ 

^(t- _ r _)( r _ r+ ) 2 (r — r_) 2 x 

2T 2 2T 2 ' 

(r - r + )(-r 2 r 2 + 4r 3 r + - 8r 2 r_r + + 6rr 2 r + - 3r 2 r 2 + 6rr_r 2 - 4r 2 r 2 ) 



Q 2 (r-r_)(r-r + ) utt 
2r 2 



8r 6 (r - r_) 



//U 2 (r — r + )(2r 2 — rr_ — 3rr+ + 2r_r + ) 0(2r — r_)(r — r + ) . tr 
+ I 40 r 3( r _ r _) + ) h 

, iQ/j,(r — r_)(r — r+) i//(r_ — r + )(r — r + )(2r 2 — rr_ — 3rr + + 2r_r + )x t2 
r 2 40r 5 (r — r_) 

(-16r 3 r_ + llr 2 r 2 + 34r 2 r_r + - 24rr 2 r + - r 2 r\ - 16rr_r 2 + Ylr 2 _r\ 
+ I 8r 4 (r-r_)(r-r+) 

/ u 2 r_(r — r + )(— 2r 2 + rr_ + 3rr + — 2r_r + ) ,u 2 2 (r — r_)x 
4fi 2 r 4 (r-r_) 2 + Y + 2(r - r+) ' 

Afi(— 4r + 3r_)(r — r+) Z/U 3 (r — r_|_)(2r 2 — rr_ — 3rr_|_ + 2r_r + ) x 
^ 2r 3 40 2 r 3 (r — r_) 

(y. - r _)( r - r+ )(-rr_ - rr + + 2r_r + ) ^ 2 (r - r_)(r - r+) , 
+ I 4 r 6 + 2r 2 ' 

2( r - r _)(-3r + 2r_)(r - r + ) 2 y?(r - r + ) 2 (2r 2 - rr_ - 3rr + + 2r_r + ) , 
+ ( ^5 + > nf 

H 2 {r - r + ) 2 (-2r 2 + rr_ + 3rr + - 2r_r + ) 
+ 80V nh 

( r — r_) 2 (r — r + ) 2 / u 2 (r — r + ) 2 (2r 2 — rr_ — 3rr + + 2r_r + )x 
+ ( ^ + 40275 )*h» = 

(5.5) 
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The initial data (as r — > oo) must obey this constraint, and then is automatically preserved 
along the integration. We have used it to verify that the integration errors were indeed 
small. 

In the previous sections we deduced the existence of an instability by finding zeros 
of a ratio, R, of the irregular to regular solutions near the horizon. As it happened, the 
functions were essentially independent near the horizon, so we could investigate their ratios 
independently. Changing the dilaton perturbation, say, and thus Rdu, did not affect the 
ratios for the other perturbations, therefore we could vary our parameters freely to find the 
instability. However, when the charge becomes large, this is no longer the case, changing 
the initial value for the dilaton perturbation will now change the ratio of field h zz , say. 
The task of finding where a 5-dimensional vector (corresponding to the ratio of each of 
the perturbation fields) is much harder. We have thus resorted to the Newton-Raphson 
method described in [17]. In a nutshell the idea is to assume that we are not too far from 
the zero 

r) /? 

= Ri(x + 6x) « Rdx) + — -Sxj (5.6) 

dxj 

We can then invert this equation to get 

ad —i 

Sx^^P- Ri{x) (5.7) 

dXj 

that is we calculate the gradient of R at the value of x we are and follow it in towards the 
zero of R. 

Figure 5 shows O for an unstable mode as a function of the charge for different values 
of \i 1 the frequency in the orthogonal direction. The important point to note is that the 
curves cross zero, at which point the unstable mode disappears for the appropriate value 
of [l. 

We have checked using the gauge invariant variables of appendix A that the solutions 
found correpsond to real physical solution and not gauge artifacts. This is obvious for the 
small charge case as the dominant quantities entering in equations (A. 4) and (A. 5) are the 
chargeless terms with the small charge terms being sub-dominant. For large charges we 
have calculated the quantity (A. 4) directly. It is non-zero for the perturbations studied. 

27 



Near the extreme caase, the function calculating the ratio between the regular and 
irregular modes has very large derivatives and thus is not easily amenable to numerical 
study. As this happens to be also at very small values of O we conjecture that the extreme 
black hole is stable. This important particular case will be investigated elsewhere [18]. 

6) Conclusion. 

We have investigated the instability of 'magnetically' charged dilatonic black p-branes 
in string theory. In particular we have shown that the chargeless instability found in [13] 
still persists in the presence of a charge on the black string or brane. This demonstrates 
that this instability is a generic phenomena. We remind the reader that this instability can 
be pictured as the string horizon collapsing in some regions of the extra dimensions and 
expanding in others. The unstable mode remains present for large values of charge but in 
this paper we have not proved its existence or otherwise for the extremal case; from the 
trend of the time frequency O as function of the charge it seems rather plausible that this 
case might be stable. However, without a more detailed analysis tailored to the extremal 
case, we cannot be definite. We can however indicate a different way of stabilizing the black 
p-brane - compactification. Compactification implies that the values of pn are quantized. If 
the compactification is on a scale smaller than the inverse mass of the black hole, the first 
allowable value of /i would be greater than the maximum one allowed for the instability, 
so such black doughnuts would be stable. Since there must be compactification of any 
extra dimensions on an extremely small scale, all but the tiniest black doughnuts would 
be safe, and those that would not would presumably have evaporated producing their own 
naked singularities long ago. Thus this instability will have no effect for contemporary 
astrophysical black holes. 

One of the main questions of interest about this instability is the nature of its endpoint 
- what does the black string become? Of course, since our calculation is linear, we cannot 
strictly say anything about the final state, however, since in the chargeless case there are no 
other scales in the problem except the mass per unit length of the black string, it seems very 
reasonable to conjecture that the endpoint of the instability will be the fragmentation of the 
black strings into a bunch of small spherical black holes. This suggestion is supported by 
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the entropy argument given in the first section. Periodic black hole solutions are known 19 , 
so unlike the Reissner- Nordstrom instability, there is a final state solution here. It is quite 
likely that such a solution will itself exhibit a Jeans-like instability to long wavelength 
clustering of black holes which itself will be unstable and so on. In any case, a process 
of fragmentation will produce a naked singularity and hence violate cosmic censorship, 
however, we should again stress that this is speculative; ideally one should follow the 
instability numerically into the non-linear regime to see if fragmentation develops, this 
work is currently underway. 

For the charged black p-brane, there are other possibilities. The fact that 5F = 
suggests that where the horizon shrinks the solution becomes locally closer to extremality. 
It is possible that when a charge is present on the black string or brane fragmentation will 
not occur. The electromagnetic repulsion might stabilize the string as the 2-spheres where 
the charge resides shrink to smaller volume. Unlike the chargeless case, the charged case 
might not lead to naked singularities and fragmentation. This conclusion assumes both 
the stability of extremal black p-branes, which is likely, as well as the non-participation of 
charge in the instability, which is an unknown without following the instability into the 
non-linear regime. Alternatively, the 'magnetic' nature of the charge considered also makes 
it plausible that as the horizon shrinks in some regions, higher energy physics might come 
into play, and that in such regions, the relevant monopole solution could 'pop out' from 
behind the horizon. For example, consider D = 4. Here an instability is known for the four 
dimensional monopole black hole, in which for small enough horizons, the t'Hooft Polyakov 
monopole is the preferred exterior solution 20 . Thus, in the case of the 6- brane should the 
horizon volume become too small, there will be sections of the rippled brane that will have 
t'Hooft Polyakov type monopole fields surrounding them. The charge, having appeared 
in an exterior smooth form, would then presumably pose no obstruction to the interior 
collapse of the event horizon, and thence a naked singularity. Again, this assumes some 
knowledge of the non-linear regime. 

We can also comment on black p-branes with topological charge. For example, consider 
the axionic black holes of Bowick et. al. 21 These carry 'quantum' charge, detectable only by 
an Aharanov-Bohm scattering process, but as far as classical physics is concerned (hence 
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our instability) they behave as if uncharged, since the quantum charge they carry does 
not affect the exterior spacetime which is consequently Schwarzschild. The field strength 
of the axion field is zero throughout the spacetime, but the gauge field (rather like the 
gauge field of a local cosmic string) is non-trivial due to the topology of the spacetime, 
B = ^gd9 A d(f). This solution can clearly be extended to a five dimensional black string 
(or ten-dimensional black 6-brane). This solution too will be unstable, however, a five (or 
ten) dimensional black hole cannot carry the same type of axion charge. Drawing analogy 
to the gauge field of a local cosmic string, it seems likely that during the fragmentation 
process, higher energy physics could come into play, producing an axion vortex, which 
could appear from behind the event horizon by an analogous process discussed above for 
the charged black p-brane. The endpoint in this case might be a line of black holes threaded 
by a cosmic axion string (not to be confused with the four-dimensional global string). 

To summarize, whether charged or uncharged, black p-branes exhibit a long wave- 
length instability which causes the apparent horizon to ripple. A plausible endpoint of this 
instability is fragmentation, and hence violation of cosmic censorship, however, without 
more evidence in the non-linear regime, this remains speculative. Perhaps a more realistic 
conclusion is that due to this instability, black strings and p-branes will not form from 
collapse in the first place. From our work, it appears there are two ways of avoiding this 
instability. One is if the black brane is extremal, although the methods described here 
were unable to probe this end of parameter space, all indications were that these solutions 
were stable. This is also anticipated since extremal branes are supersymmetric. The other 
way of avoiding the instability is by compactification, indeed, to be totally speculative, 
since the endpoint of the instability does appear to be some periodic solution, it is tempt- 
ing to suggest that the instability of some primordial black string could have triggered an 
effective compactification of our universe! 

Acknowledgments. 

We are grateful to J.B.Hartle, J.A.Harvey, S.W. Hawking, G.Horowitz, D. Wiltshire, 
and W.H.Zurek for useful conversations. We appreciate the careful perusal of the manuscript 
and accompanying remarks by J.Bowcock and P.Laflamme. R.G. is supported by the 

30 



S.E.R.C. R.L. thanks Los Alamos National Laboratory for support. R.G. also acknowl- 
edges the support of the McCormick Fellowship at the Enrico Fermi Institute where this 
work was started. 

References. 

[1] S.W.Hawking, Commun. Math. Phys. 43 199 (1975). 

[2] S.W.Hawking, Commun. Math. Phys. 87 395 (1982). 

[3] E.Poisson and W. Israel, Phys. Rev. D41 1796 (1990). 

[4] G.Gibbons and K.Maeda, Nucl. Phys. B298 741 (1988). 

D.Garfmkle, G.Horowitz and A.Strominger, Phys. Rev. D43 3140 (1991). 

[5] R.Gregory and J.Harvey, Phys. Rev. D47 2411 (1993). 
J.Horne and G.Horowitz, Nucl. Phys. B399 169 (1993). 

[6] R.Myers and M.J.Perry, Ann. Phys. (N. Y.) 172 304 (1986). 

[7] G.T.Horowitz and A.Strominger, Nucl. Phys. B360 197 (1991). 

[8] W.Israel, Commun. Math. Phys. 8 245 (1968). 

S.W.Hawking, Commun. Math. Phys. 25 152 (1972). 

[9] R.Gregory and R.Laflamme, Phys. Rev. D37 305 (1988). 
[10] C.V. Vishveshwara, Phys. Rev. Dl 2870 (1970). 
[11] P.Bizon and RWald, Phys. Lett. 267B 173 (1991). 
[12] B. Whitt, Ph.D Thesis, Cambridge, 1988. 
[13] R.Gregory and R.Laflamme, Phys. Rev. Lett. 70 2837 (1993). 

[14] S.W.Hawking and G. F.R.Ellis, The Large Scale Structure of Spacetime (Cambridge 

University Press 1973). 
[15] R.M.Wald General Relativity (Chicago University Press). 
[16] J.Bardeen, Phys. Rev. D22 1882 (1980). 

[17] W.H. Press, B.P.Flemming, S.A.Teukolsky and W.T. Vettering, Numerical Receives 

(Cambridge University Press 1988). 
[18] R.Gregory and R.Laflamme, The stability analysis of the extremally charged black 

string, in preparation. 
[19] A.Bogojevic and L.Perivolaropoulos, Mod. Phys. Lett. A6 369 (1991). 

31 



[20] K.Lee and E.Weinberg, Phys. Rev. Lett. 68 1100 (1992). 

[21] M.Bowick, S.Giddings, J.Harvey, G.Horowitz and A.Strominger, Phys. Rev. Lett. 61 
2823 (1988). 

Appendix A. Gauge invariant variables. 

One of the worries of results in linear perturbations of a system with gauge invariance 
is to insure that the result are not artifact of the particular gauge chosen. One easy way 
to deal with this problem is to use gauge invariant variables. In this appendix we explain 
how to construct these variables. 

Under a coordinate transformation 

x a ->x a +£ a (A.l) 
we have the following transformation of the fields 

h ab _^ h ab + £(a;6) 

0^0 + £ a L> a (A2) 

fiF'ai...a D -2 = Fa 1 ...a D _ 2 ;d£, + (D — 2)i^[o2...0D_ 2 ^;oi] 

In the metric given in eq. (2.3-2.5), they have the explicit form (for V+ = (r — r+) and 

V- = (r — r_)) 



t V V+ , V 
h rt - h rt + -^±?' - ^ne r (A3c) 



fc** - fc** - + ^ (A3e) 

fc" - /i 22 + (A.3g) 
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It is easy to show that the linear combination 

Y = h* r + i^rh zz ' - ^h" (AA) 
2\xr l 2 

is indeed invariant under gauge transformation. However this gauge invariant variable is 
not well defined when // = 0, a case under consideration in section 3. A more complicated 
variable 

x h tr | (r+-R-)V+ tt VI tt , 

, (A5) 

(r+ — r_)(— r + 2rr+ + 2rr_ — 3r+r_) (r+ — R-)r , V-flr 

— K H — K K 

V+V-Q 40 2V+ 

can be used to investigate the \x — > limit in the chargeless case. In this limit we get 
I « ft as can be shown from the asymptotic region (r — > oo and r — > r + ) discussed in 
section 3. This analysis shows that the perturbations described in sections 3-5 describe 
real physical instabilities except in the case fx = 0. 

Appendix B) Background Quantities. 

The background metric is given explicitly by the following: 

rD _3_ r D-3 
9tt ~ r D-3 _ r ^"3 

r 2(£>-3) 

9ri 



(j.D-3 _ r D-3^ r D-3 _ r ^~ 3> 

a 

9e a e =r 2 S a p Y\ sin2 °n for a > 1 

n=l 



(BA) 



9ij = S 



t-3 



where a, f3 run from 1 to D — 2, and represent the angular variables. The 'matter' content 
has solution: 

•=-H i -(Tn (b - 2) 
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F = Qeo-2 , Q 2 = ^r^(r + r_) D - 3 



The form of the perturbation is 



g ac g bd 5g cd = h ab = e nt+ ^ ixl 



H ij (r) 


H it 


H ir 







H tt 


R tr 







H tr 


jijrr 














K(r) 




















Kj sin 2 6 



(B.3) 



(BA) 



and finally 



5$ = e nt+iliixi f(r) 



5F = 

The non-zero elements of the Riemann tensor are given by: 



(B.5) 



(B.6) 



Bj 



rtr 



JDS) 


"(^) D - 3 - C-) D ~ 3 


~D-2-C^r) D -* 


2r 2 






^-( r -) D ~\ 


2 



R 



9 a lphatd p 



n r9 fj r 



2r 2 



r r 



9otf3 



_-(£>- 3) 


(- 


) D - 3 + (^) D 


-3 


-2(^)°- 3 ] 


2r 2 


1-(^) D - 3 " 




1 





Ricci tensor 



(S.7a) 

(B.76) 
(5.7c) 
(S.7d) 



-Rtt — 



(D-3) 2 /r_\^"3 



2r 2 



r+ \D-3 



l-(^) 



D-3 



(S.8a) 



(-D - 3)(^)^" 3 [(3D - 7)(^)°- 3 + (D - l)( r -) D - s - 2(D - 2) (l + ( r -^) D ~ 3 )] 



2r 2 [l-(I-) 
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(D-3) 



(t) 



D-3 



9a/3 



(5.86) 
(5.8c) 



Ricci scalar 



* = ^ a -i^r^ i-(g-^-(g^)(v) p - 3 (v» o " 3+ ' 2£, - 5 '' ! 7 )p " 3 - (! 7 )g " 31 

(B.9) 

The gauge conditions component by component are given by the transverse conditions 
V ' a h ah = 0. For simplicity in what follows, we will write R± = (r±/r) D ~ 3 

+ (D -yp-^ + ^g^ (r) _ 

„ - h ar _ (D - - R+)(R+ - R-) htt _ [(D - 2) - (25 - 5)5_] 
Va ^ " r(l-5_) * r(l-5_) * 

Cl — i? I 2 

- 2(5 - 2)r(l - 5_)(1 - 5 + )if + fi/i tr + i A t i /i r< + ^ ±^h M ' 

V a /i ae = V a /i a0 = (5.10c) 

v hai = m(l-R + ) h u + nh u _ m h rr + [(5-2) - (25-5)5_] 

2(1 -5_) 2(l-5 + )(l-5_) r(l-5_) 

_ (5-2) .^2 K + ^ _ + ^ ri / (r) = Q 

2 2 7 

(5.10d) 



Finally we give equations for the perturbations /i ab , the trace /i = g a bh ab and / = 5<E>. 

-(5 - 3)/i rr (5+ - R-)[2(D - 2) - (5 - 1)5+ - (35 - 7)5_ + 2(5 - 2)5+5_] 
" ^2 (l-5+) 3 (l-5_) 
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, h «( f|2 n 2 (i-R.) (D-3) 2 (R+-R-) 2 \. m 9Vn . (R+-R-)(i-R-) 

2(D - 3)0 (R+ + R+R- - 2R-) ^ ^ 2 (1 - i?-) 2 /(r) | 2(£> - 3) (1 - R-)(R+ - g_) ;/ 



r (l-^+) 2 (l- J R + ) 2 ' /w r 

lit' 

+ [(L>-2)-(2 J D-5) J R_ + (2L>-7) J R + -( J D-4) J R +J R_] + (l- J R + )(l-i?_)/i" (S.llo) 

r 



2 2 (1- J R_) 4( J D-2)(l-2i? + -(D-l) J R_) 



+ 



(1-12+) 2r*(l- J R+)(l- J R_) 

(3D 2 - 26D + 47)R 2 _ + (D 2 - 10 D + 17)R 2 + + 2(2D 2 - IAD + 23)R 2 + R 2 _ 

2r 2 (l- R+)(l- R-) 

2(L> 2 - 5)R + R_ + (D 2 - 12D + 23)R 2 + R_ + (5D - 13) (D - 5)R 2 _R+ \ 

r 2 (l - R+)(l - R-) ) 

+ (D - 2)K(1 - R+)(l - R_)[2 + 2(D - 2)R+R_ - (D - 1)(R+ + R_)] 
2(D-3)Sl{R+-R-) htr + [(D _ 2) _ _ 5)R+ _ (4jD _ U)R _ + (5jD _ u)R + R_]^— 



2( - D 3 \ l-R + )(l-R-)[R + +R--2R + R.]f'-4(l-R + ) 2 (l-R.) 2 f"(r)+(l-R+)(l-R-)h rr " 



(B.llb) 

>h rt 

<KD ~ 3) f hrr [2R--R + -R + R-] _ h « _ \ _ mf (D - S)(R + - fl_)(l - R_) 



(1-R+) 2 v ^ -'J J r(l-R+) 



+ h rt | _f _ ^ 2 (1 " R-) (D - 2)[1 - 2(D - 2)tf_ + (D — 5)R+] 



(1-R+) r 2 (l-R + )(l-R-) 

(D 2 - 10D + 20)R 2 _ -R\- {2D - 5))R 2 + R 2 _ 
r 2 (l- R+)(l- R-) 
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(3D 2 - 17 D + 26)R 2 + R_ - (3D 2 - 13D + 16)R + R_ - (2D 2 - 18D + 34)R 2 _R. 
+ r 2 (l -#+)(!-#_) 



hrt' 

+4n(l-R_) 2 f' + l(D-2)-R + -(3D-8)R_ + (2D-5)R + R-] + (1 -R+)(l -R-)h r 



(B.llc) 

•h ee 

_ h rr [2 — (D — 1)R+ - (3D - 7)R- + 2(2D - 5)R+R-] (D - 3)h tt (R+ - R-)(l - R+) 
7^ (1-R+)(1-R.) + 7 A (i-R-) 

+K[ V - + - 2) + (D - 4)(R + + R_) - (D 2 -3D- 1)R + R-}) 

-4, f\r) (1 - - +[(D+2)-5(R + +R.)-(D-8)R + R.]^ + (l-R+)(l-R.)K" 

(B.lld) 

•h tl 

„ _ (D - 3M R + - 2H- + fl + B-) h „ _ m(D-3)R tr + (1 - H-) 



r(l- J R + ) 2 r(l-i?_) ^ JW (1- J R + ) 



fr*^ "^ 1 (Rile) 



(l-i2+) V r 



•h ri 



= _( J -3)ncH + -fl-) h , < _ WJ-3)«- r _ W(r)(1 _ fl+)(1 _ B _, 

2 (1 - i?_) (1 - R+)[-(D - 2)(1 — (-D — l)i?_) - (2D - 5)# 2 ]\ 



(1-/2+) r 2 (l-i?_) / 

(1 -R+)(D-2- (2D - 5)R_) + (1 - R+)(l - R_)h ri " (5.11/) 
r 

(^-3)«-r,.. ,H , .... , ... .. „ x ,»Y..2 , fi 2 (l-^-)' 



r(l - R ) [i » ihrJ + + A ^ f{r) " ^ + (1-R + ) J 



= 

r(l - R-) L ' ~ ' " J 
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+ — [(D — 2) — R- — R + — (D — 4)R + R_] + (1 - R+)(l - R_)h 13 (B.llg) 
For the trace, h, 

= -2(D -2){D- 3) 2 R + R-K - (> + ^rffi ) ( h ~ 4 ^ 

+ (l-R-)(D-2-R+) _ , _ _ „ _ „ (B.Uh) 

r 

and finally for the dilaton perturbation: 

(D -3) 2 R_(1- R + )(R + -R-) ltf (D-3)(D-2) „ n . /r% , „ . 
= 4^(1-^)2 + 2 KR - [1 ~ {D ~ 2)R+] 

(D - 3)R_[2{D - 2) - (3D - 7)(R+ + R_) + 2(2D - 5)R + R_] / 2 Q 2 (l - 

4r 2 (f-#_) 2 (f-#+) (1-R+) ) 

f 

+ — [(D — 2) — R + + (D — A)R_ - {2D - 7)R + R_] + (f - R+)(l - R_)f" (B.lli) 
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Captions. 

Figure 1. Penrose diagram for the metric of eq.(l). All but two dimensions have been 
suppressed; each point of this diagram corresponds to a D-2 sphere times 10-D flat space. 
The line at v = vq is the initial hypersurface on which the perturbations are set. It is a 
Cauchy surface for the spacetime exterior of the black hole. 

Figure 2. Plot of the modes H u , H tr and H rr as a function of (r — r+) for a regular mode. 
As (r — r_|_) becomes large all the functions decay exponentially, they reach a maximum at 
small (r — r + ) and go to zero as (r — r + ) — > 0. 

Figure 3. Plot of O as a function of \± for black strings and branes with D = 4, ...,9 
and r_|_ = 2 for which an instability has been found. The bold points correspond to value 
calculated numerically and the lines have been traced to guide the eye. 

Figure 4. Schematic diagram of the instability. All the dimensions have been suppressed 
except the r and z one. The geometry is initially invariant under a z translation. The 
instability increases the size of the apparent horizon at some values of z and decreases it 
at other values. In the chargeless case it seems reasonable that the apparent horizon will 
brake into different parts. 

Figure 5. Plot of the modes /, q, H zz and H tz as a function of (r — r+) for a regular mode. 
As (r — r+) becomes large all the functions decay exponentially, they reach a maximum at 
small (r — r + ) and go to zero as (r — r + ) — > 0. 

Figure 6. Plot of O as a function of [i for a charged 5d black string with charges corre- 
sponding to r_ = 0, 1.0, 1.5 and r + = 2 for which an instability has been found. The bold 
points correspond to value calculated numerically and the lines have been traced to guide 
the eye. 
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